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Abstract 

In this paper, we study the SIS (susceptible-infected-susceptible) and SIR (susceptible- 
infected-removed) epidemic models on undirected, weighted networks by deriving pairwise- 
type approximate models coupled with individual-based network simulation. Two dif- 
ferent types of theoretical/synthetic weighted network models are considered. Both 
models start from non-weighted networks with fixed topology followed by the alloca- 
tion of link weights in either (i) random or (ii) fixed/deterministic way. The pairwise 
models are formulated for a general discrete distribution of weights, and these models 
are then used in conjunction with network simulation to evaluate the impact of dif- 
ferent weight distributions on epidemic threshold and dynamics in general. For the 
SIR dynamics, the basic reproductive ratio i?o is computed, and we show that (i) for 
both network models Rq is maximised if all weights are equal, and (ii) when the two 
models are equally matched, the networks with a random weight distribution give rise 
to a higher Rq value. The models are also used to explore the agreement between the 
pairwise and simulation models for different parameter combinations. 

Keywords: weighted network, pairwise model, moment closure 
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1 Introduction 

Conventional models of epidemic spread consider a host population of identical individuals, 
each interacting in the same way with each of the others (see [2j [2H H5] and references 
therein). At the same time, in order to develop more realistic mathematical models for 
the spread of infectious diseases, it is important to obtain the best possible representation 
of the corresponding transmission mechanism. To achieve this, more recent models have 
included some of the many complexities that have been observed in mixing patterns. One 
such approach consists in splitting the population into a set of different subgroups, each 
with different social behaviours. Even more detail is included within network approaches, 
which allow to include differences between individuals, not just between sub-populations. 
In such models, each individual is represented as a node, and interactions that could permit 
the transmission of infection appear as edges linking nodes. The last decade has seen a 
substantial increase in the research of how infectious diseases are spread over large networks 
of connected nodes, where networks themselves can represent either small social contact 
networks or larger scale travel networks, including global aviation networks [22| |26| \37\ H6| 
1521 1531 EH1 ED]- Importantly, the characteristics of the network, such as the average 
degree and the node degree distribution have a profound effect on the dynamics of the 
infectious disease spread, and hence significant efforts are made to capture properties of 
realistic contact networks. 

Network models provide an intuitively appealing way to consider the social structure of 
populations, but they bring with them the challenge of collecting sufficient data to parame- 
terise them fully. Ideally, an entire population would be sampled and all relevant interactions 
measured in order to reconstruct the contact network, but this is, unsurprisingly, seldom 
possible, and it is generally necessary to make suitable approximations. With the advance 
of computer-based tracking, statistical properties of many realistic networks (e.g. mobile 
phone calls) have been studied in great detail, and these are often used as proxy for the 
analysis of epidemic dynamics ^2E\ E] . 

One of the simplifying assumptions often put into models on networks is that all links are 
equally likely to transmit infection [TH |3TJ H6j [63] . However, a more detailed consideration 
leads to an observation that this is often not the case, as some links are likely to be far more 
capable of transmitting infection than others due to closer contacts (e.g. within households 
[TT] ) or long-duration interactions [JSJ USB E2J E2J IM] . To account for this heterogeneity in 
properties of social interactions, network models can be adapted, thus resulting in weighted 
contact networks, where connections between different nodes have different weights. These 
weights may be associated with the duration, proximity, or social setting of the interac- 
tion, and the key point is that they are expected to be correlated with the risk of disease 
transmission. The precise relationship between the properties of an interaction and its risk- 
iness is hugely complex; here, we will consider a "weight" that is exactly proportional to the 
transmission rate along a link. From a general perspective, one can consider two types of 
weighted networks: symmetric networks, where the strength only varies between different 
edges but is the same for a given pair of connected nodes, or go one step further and allow 
for the strength to be different even for the same pair of connected nodes A and B, depend- 
ing on whether A infects B, or B infects A. The first type of networks can be efficiently 
conceived as undirected graphs represented by symmetric connection matrices, while for the 
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second one it is necessary to consider the contact network as a directed graph, resulting in a 
connection matrix, which in general is non-symmetric. Although consideration of weighted 
networks may seem as an additional complication for the analysis of epidemic dynamics, in 
fact it provides a much more realistic representation of actual contact networks. 

Early models of weighted networks arose in the context of identifying genetic and metabolic 
networks [TJ EE] . Later on, they found applications in a number of diverse research areas, 
such as international trade networks [TQJ |32], scientific collaboration networks [HI E21 EI], 
aviation networks of individual countries [3J EH] or global aviation networks [5J EDI SI] and 
networks of mobile phone communications [58]. Yook et al. [73 J started a systematic study 
of evolving networks with nonbinary connectivities by considering growing weighted and 
exponential networks. This approach was later developed into a generic formalism for the 
analysis of evolving weighted networks [51 El El EDI EI] ■ Newman [55] has discussed general 
properties of weighted networks. Substantial amount of work has been done on the analysis 
of scale-free networks with different types of weight distribution [671 [73]. Wang et al. [68] 
have proposed a model with dynamical adjustment of weights in technological networks. 

In epidemiological context, Yan et al. [70] analysed a model on weighted scale-free 
networks and found that heterogeneity in weight distribution leads to a slowdown in the 
spread of epidemics. Furthermore, they have shown that for a given network topology and 
mean infectivity, epidemics spread fastest in unweighted networks. Chu et al. [17] have 
investigated the dynamics of an SI model on weighted scale-free networks with community 
structure and showed that the contribution coming from weights of edges connecting different 
communities exceeds the impact of weights on edges inside communities. Karsai et al. [33] 
have studied mean-field and finite-size scaling in weighted networks with small weights on 
edges connecting high-degree nodes. Yang et al. [72] have shown that disease prevalence can 
be maximized when the edge weights are chosen to be inversely proportional to the degrees 
of receiving nodes but, in this case, the transmissibility was not directly proportional to the 
weights and weights were also asymmetric. Yang & Zhou [71] have considered SIS epidemics 
on homogeneous networks with uniform or power-law edge weight distribution and shown 
how to derive a certain type of mean-field description for such models. Britton et al. [T5] 
have derived an expression for the basic reproductive ratio in weighted networks with generic 
distributions of node degree and weights, and Deijfen [23] has performed a similar analysis 
to study vaccination in such networks. In terms of practical epidemiological applications, 
weighted networks have already been effectively used to study control of global pandemics 
[HI Ell EH] and the spread of animal disease due to cattle movement between farms [351 EH] ■ 
Eames et al. [29] have considered an SIR model on an undirected weighted network, where 
rather than using some theoretical formalism to generate an idealized network, the authors 
have relied on social mixing data obtained from questionnaires completed by members of 
a peer group [62] to construct a realistic weighted network. Having analysed the dynamics 
of epidemic spread in a such a network, they showed how information about node-specific 
infection risk can be used to develop targeted preventative vaccination strategies. As an 
alternative approach to modelling heterogeneity in network interactions, Eames [27] has 
considered a network with two types of contacts: random and regular, where the first 
type refers to arbitrary connections with any node in the network, while the second type 
designates multiple repeated contacts with the same nodes. He showed that in a highly 
clustered population, random contacts allow infection to reach otherwise inaccessible parts 
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of the network, while in the case of all contacts being regular, clustering leads to a significant 
reduction in the spread of infection. 

In this paper, we consider the dynamics of an infectious disease spreading on weighted 
networks with different weight distributions. Since we are primarily concerned with the 
effects of weight distribution on the disease dynamics, the connection matrix will be assumed 
to be symmetric, representing the situation when the weights can only be different for 
different network edges, but for a given edge the weight is the same irrespective of the 
direction of infection. From epidemiological perspective, we consider both the case when 
the disease confers permanent immunity (represented by an SIR model), and the case 
when the immunity is short-lived, and upon recovery the individuals return to the class 
of susceptibles (SIS model). For both of these cases we derive the corresponding ODE- 
based pairwise models and their closure approximations. Numerical simulation of both the 
network dynamics and the pairwise approximations are performed. 

The outline of this paper is as follows. In the next section, the construction of specific 
weighted networks to be used for analysis of epidemic dynamics is discussed. This is comple- 
mented by the derivation of corresponding pairwise models and their closure approximations. 
Section 3 contains the derivation of the basic reproductive ratio Rq for the SIR model and 
for different weight distributions as well as numerical simulation of both network models 
and their pairwise ODE counterparts. The paper concludes in Section 4 with discussion of 
results and possible further extensions of this work. 

2 Material and methods 
2.1 Network construction 

There are two conceptually different approaches to constructing weighted networks for mod- 
elling infectious disease spread. In the first approach, there is a seed or a primitive motif, 
and the network is then grown or evolved from this initial seed according to some specific 
rules. In this method, the topology of the network is co-evolving with the distribution of 
weights on the edges [5J [7J El EH E2]- Another approach is to consider a weighted network 
as a superposition of an unweighted network with a distribution of weights across edges 
which could independent of the original network or it may be correlated with node metrics, 
such as their degree, [151 E21 EU EE]. In this paper we use the second approach in order to 
investigate the particular role played by the distribution of weights across edges, rather than 
network topology, in the dynamics of epidemic spread. Besides computational efficiency, this 
will allow us to make some analytical headway in deriving and analysing low-dimensional 
pairwise models which are likely to perform better when weights are attached according to 
the scenarios described above. 

Here we consider two different methods of assigning weights to network links: a network 
in which weights are assigned to links at random, and a network in which each node has the 
same distribution of weighted links connected to it. In reality, there is likely to be a great 
deal more structure to interaction weights, but in the absence of precise data and also for the 
purposes of developing models that allow one to explore a number of different assumptions, 
we make these simplifying approximations. 
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2.1.1 Random weight distribution 

First we consider a simple model of an undirected weighted network with N nodes where 
the weights of the links can take values Wi with probability p i: where % = 1, 2, . . . , M. The 
underlying degree distribution of the corresponding unweighted network can be chosen to be 
of the more basic forms, e.g. homogeneous random or Erdos-Renyi-type random networks. 

The generation of such networks is straightforward, and weights can be assigned during 
link creation in the unweighted network. For example, upon using the configuration model 
for generating unweighted networks, each new link will have a weight assigned to it based 
on the chosen weight distribution. This means that in a homogeneous random network 
with each node having k links, the distribution of link weights of different type will be 
multinomial, and it is given by 

k\ 

P(n wi , n W2 , . . . , n WM ) = j j ^P\ 1 P 2 2 ■ ■ -Pm > (-0 

n Wl ■Ti W2 .... ti Wm . 

where, n Wl + n W2 + • • • + n WM = k and P(n Wl ,n W2 , . . . , n WM ) stands for the probability of a 
node having n Wl , n W2 , . . . , n WM links with weights w±, W2, ■ ■ ■ , %, respectively. While the 
above expression is applicable in the most general set-up, it is worth considering the case of 
weights of only two types, where the distribution of link weights for a homogenous random 
network becomes binomial 

P(n Wl ,n W2 = k-n Wl ) =( k )p?{l - Pl ) k ~ n \ (2) 

where, pi + P2 = 1 and n Wl + n W2 = k. The average link weight in the model above can be 
easily found as 

M 

w ran d om = J2pm, 
i=l 

which for the case of weights of two types w± and w 2 reduces to 

W av r) = Pl W l + P2W2 = PlWl + (1 - Pl)w 2 . 

2.1.2 Fixed deterministic weight distribution 

As a second example we consider a network, in which each node has ki links with weight 
Wi (i = 1,2,..., M), where k\ + k 2 + • • • + k M = k. The different weights here could 
be interpreted as being associated with different types of social interaction: e.g. home, 
workplace, and leisure contacts, or physical and non-physical interactions. In this model all 
individuals are identical in terms of their connections, not only having the same number of 
links (as in the model above) but also having the same set of weights. The average weight 
in such a model is given by 

M h 
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where is the fraction of links of type i for each node. In the case of links of two types 
with weights W\ and the average weight becomes 



W aJ ] = Pl W l + V2W2 = -£Wi + -j^W 2 = yWi H 7~ w 2 



k\ ki ki k — ki 



2.1.3 Epidemic models 

In this study the simple SIS and SIR epidemic models are considered. The epidemic dy- 
namics is specified in terms of infection and recovery events. The rate of transmission across 
an unweighted edge between an infected and susceptible individual is denoted by r. This will 
then be adjusted by the weight of the link which is assumed to be directly proportional to 
the strength of the transmission along that link. Infected individuals recover independently 
of each other at rate 7. The simulation is implemented using Gillespie algorithm [SB] with 
inter-event times distributed exponentially with a rate given by the total rate of change in 
the network, with the single event to be implemented at each step being chosen at random 
and proportionally to its rate. All simulations start with most nodes being susceptible and 
with a few infected nodes chosen at random. 



2.2 Pairwise equations 

In this section we extend the classic pairwise model for unweighted networks [SI 161] to the 
case of weighted graphs with M different link-weight types. Pairwise models successfully 
interpolate between classic compartmental ODE models and full individual-based network 
simulation with the added advantage of high transparency and a good degree of analytical 
tractability. These qualities makes them an ideal tool for studying dynamical processes 
on networks [27J EH SHI HI], and they can be used on their own and/or in parallel with 
simulation. The original versions of the pairwise models have been successfully extended to 
networks with heterogenous degree distribution [28], asymmetric networks [65] and situations 
where transmission happens across different /combined routes [27, 38J as well as when taking 
into consideration network motifs of higher order than pairs and triangles [39J . The extension 
that we propose is based on the previously established precise counting procedure at the 
level of individuals, pairs and triples, as well as on a careful and systematic account of all 
possible transitions needed to derive the full set of evolution equations for singles and pairs. 
These obviously involve the precise dependency of lower order moments on higher order 
ones, e.g. the rate of change of the expected number of susceptible nodes is proportional 
to the expected number of links between a susceptible and infected node. We extend the 
previously well-established notation [H] to account for the added level of complexity due to 
different link weights. In line with this, the number of singles remains unchanged, with [A] 
denoting the number of nodes across the whole network in state A. Pairs of type A — B, 
[AB], are now broken down depending on link weights, i.e. [AB]i represents the number of 
links of type A — B with the link having weight Wi, where as before i = 1, 2, . . . , M and 
A,Be {S, I, R} if an SIR dynamics is used. As before, links are doubly counted (e.g. in 
both directions) and thus the following relations hold: [AB] m = [-BA] m and [AA] m is equal 
to twice the number of uniquely counted links of weight w m with nodes at both ends in state 
A. From this extension it follows that ~YldLi[AB]i = [AB]. The same convention holds at 
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the level of triples where [ABC] mn stands for the expected number of triples where a node in 
state B connects a node in state A and C via links of weight w m and w n , respectively. The 
weight of the link impacts on the rate of transmission across that link, and this is achieved 
by using a link-specific transmission rate equal to TWi, where i = 1,2, . . . , M. In line with 
the above, we construct two pairwise models, one for SIS and one for SIR dynamics. 
The pairwise model for the SIS dynamics can be written in the form: 

[S]=l[I]-rE A n Ll^n[SI] n , 
[i\=TY% sl W n [Sl\ n ->y[I\, 

[SI] m = 7([i7] m - [SI]m) + T Y.n=l W n ([SSI] mn - [ISI) nm ) ~ TW m [SI] m , (3) 

[77] m = -2j[II] m + 2rJ2n=i w n [ISI] nm + 2rw m [SI] m , 

[SS) m = 2j[SI) m - 2Tj2n=l W n[SSI) mn , 

where m — 1, 2, 3, M and [AB] m denotes the expected number of links with weight w m 
connecting two nodes of type A and B, respectively (A, B G {S, I}). 

In the case when upon infection individuals recover at rate 7 and once recovered they 
maintain a life-long immunity, we have the following system of equations describing the 
dynamics of a pairwise SIR model: 

[S] = -rEn=^n[SI]n, 
[I]=rEn=^n[SI} n -l[I], 

[SS] m = ~2rJ2n=l w n[SSI]mn, 

[SI) m = rJ2n=i w n ([SSI] mn - [ISI) nm ) - rw m [SI) m - j[SI) m , (4) 

[SR] m = -TY^xWnilSR^nm + j[SI) m , 
[il}m = 2rJ2 n =l Wn[ISI} nm + 2TW m [SI} m - 2j[II] m , 
[IR} m = T^nLl Wn[ISR] nm + 7 ([Ji] m - [IR] m ), 
[RR]m = l[IR]m, 

where again m — 1, 2, 3, M with the same notation as above. 

The above systems of equations ([3]) and (TJ]) are not closed, as equations for the pairs 
require knowledge of triples, and thus, equations for triples are needed. This dependency on 
higher-order moments can be curtailed by closing the equations via approximating triples in 
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terms of singles and pairs [H]. For both systems, the agreement with simulation will heavily 
depend on the precise distribution of weights across the links, the network topology, and 
the type of closures that will be used to capture essential features of network structure and 
the weight distribution. As a check and a reference back to previous model, in Appendix A 
shows how systems ^ and (jl]) reduce to the standard unweighted pairwise SIS and SIR 
models [H] when all weights are equal to each other, w\ = w 2 = ■ ■ ■ = w m = W. 



2.2.1 Closure relations 

The most natural extension of the classic closure is given by 

[ABCU = ^l\Vk, (5) 

where k is the number of links per node for a homogeneous natwork or the average nodal 
degree for networks with other than homogenous degree distributions. Even for the simplest 
case of homogenous random networks with two weights, a different version of the closure can 
be derived. The starting point for this different closure is the observation that the average 
number of links of weight w\ across the whole network is k\ = pik < k, and similarly, the 
average number of links of weight w 2 is k 2 = (1 — Pi)k < k. This motivates the following 
closures 

[BC}i h-liAB^BCl 



[ABC] ii = [ABUh - 1 
[ABC] 12 = \AB]ik 2 
\ABC] 2l = [AB] 2 k x 



[ABC} 22 = [AB} 2 (k 2 - 1) 



h[B] h [B] 

[BC} 2 _ [ABUBC], 
k 2 [B] - [B] ' 

[BC]i _ [ABUBC], 
h[B] - [B] ' 

[BC} 2 k 2 - 1 [AB\ 2 [BC\, 



(6) 



k 2 [B\ k 2 [B\ 

The specific choice of closure will depend on the structure of the network and, especially, how 
the weights are distributed. For example, for the case of the homogeneous random networks 
with links allocate randomly, both closures offer a viable alternative. For the case of a 
network where each node has a fixed pre-allocated number of links with different weights, 
e.g. ki and k 2 links with weights Wi and w 2 , respectively, the second closure §6§ offers the 
more natural/intuitive avenue towards closing the system and obtaining good agreement 
with network simulation. 



3 Results 

In this section we present analytical and numerical results for weighted networks and pairwise 
representations of SIS and SIR models in the case of two different link-weight types (i.e. 
Wi and u^)- 
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3.1 Threshold dynamics for the SIR model - the network per- 
spective 

The basic reproductive ratio, R (the average number of secondary cases produced by a 
typical index case in an otherwise susceptible population), is one of the most fundamental 
quantities in epidemiology ([21 [25]). Besides informing us on whether a particular disease 
will spread in a population, as well as quantifying the severity of an epidemic outbreak, it 
can be also used to calculate a number of other important quantities that have good intuitive 
interpretation. In what follows, we will compute Rq and Ro-\ike quantities and will discuss 
their relation to each other, and also issues around these being model-dependent. First, we 
compute Rq from an individual-based or network perspective by employing the next gener- 
ation matrix approach as used in the context of models with multiple transmission routes 
such as household models 



Random weight distribution: First we derive an expression for R when the underlying 
network is homogeneous, and the weights of the links are assigned at random according to a 
prescribed weight distribution. In the spirit of the proposed approach, the next generation 
matrix can be easily computed to yield 

(k — V)P\Ti {k — l)pi7"i 

(k - l)p 2 r 2 (k - l)p 2 r 2 



NGM = ( aij ) iij=1 , 2 

where 



TWi TW 2 

n = ■ — , r 2 



TW\ + 7 TW 2 + 7 

represent the probability of transmission from an infected to a susceptible across a link 
of weight ui\ and w 2 , respectively. Here, the entry stands for the average number of 
infections produced via links of type i (i.e. with weight Wi) by a typical infectious node who 
itself has been infected across a link of type j (i.e. with weight Wj). Using the fact that 
p 2 = 1 — pi, the basic reproductive ratio can be found from the leading eigenvalue of the 
NGM matrix as follows 

i2S = (A;-l)(p 1 ri + (l-pi)r 2 ). (7) 
In fact, the expression for Rq can be simply generalised to more than two weights to give 
Ro = (k — 1) J2tLi P% r %-, where w m has frequency given by p m with the constraint that 
J2i=iPi = 1- ft i s straightforward to show that upon assuming uniform weight distribution 
Wi = W for i = 1, 2, . . . , M, the basic reproduction number on a homogeneous graph reduces 
to Rq = (k — l)r as expected, and where, r = tW/(tW + 7). 

Deterministic weight distribution: The case when the number of links with given weights for 
each node is fixed can be captured with the same approach, and the next generation matrix 
can be constructed as follows 

NGM = ( kl ~ 1)ri kiri 

k 2 r 2 (k 2 - l)r 2 

As before, the leading eigenvalue of the NGM matrix yields the basic reproductive ratio, 



d2 (k, - l)n + (k 2 - l)r 2 + y/lfa - Vyr\ - (k 2 - l)r 2 ] 2 + ^k x k 2 r x r 2 

K - o ■ 
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Using these two equivalent expressions for the basic reproductive ratio, it is possible to prove 
the following result. 

Theorem 1. Given the setup for the fixed weight distribution and using p± = k\jk, 
P2 = kijk and kx + ki = k, if 1 < k\ < k — 1 (which implies that 1 < k2 < k — 1), 
then Rl<R l Q . 

The proof of this result is sketched out in Appendix B. This Theorem effectively states that 
provided each node has at least one link of type 1 and one link of type 2, then independently 
of disease parameters, it follows that the basic reproductive ratio as computed from ([7]) 
always exceeds or is equal to an equivalent R computed from ([8]). 
It is worth noting that both Rq values reduce to 



Rl = Rl = R n = (k-l)r= {k J )TW . 



^--o-Ro = (k-l)r = ^ w f T , (9) 

if one assumes that weights are equal, i.e. ui\ = wi = W. As one would expect, the first 
good indicator of the impact of weights on the epidemic dynamics will be the average weight. 
Hence, it is worth considering the problem of maximising the values i?o under assumption 
of a fixed average weight: 

PiWi + p 2 w 2 = W. (10) 
Under this constraint the following statement holds. 

Theorem 2. For weights constrained by p\W\ + P2W2 = W (or (ki/k)w\ + (k2/k)w2 = W 

for a fixed weights distribution), R$ and Rq attain their maxima when w\ = W2 = W , and 

(h _ iWvi/ 

the maximum values for both is Rn = (k — l)r = — — . 

v ' rPU + 7 

The proof of this result is presented in Appendix C. 

The above results suggest that for the same average link weight and when the one-to- 
one correspondence between p\ and ki/k, and P2 and /^/^ holds, the basic reproductive 
ratio is higher on networks with random weight distribution than on networks with a fixed 
weight distribution. This, however, does not preclude the possibility of having a network 
with random weight distribution with smaller average weight exhibiting an i?o value that 
it is bigger than the Rq value corresponding to a network where weights are fixed and the 
average weight is higher. The direct implication is that it is not sufficient to know just 
the average link weight in order to draw conclusions about possible epidemic outbreaks on 
weighted networks; rather one has to know the precise weight distribution that provides a 
given average weight. 

Figure [T] shows how the basic reproductive ratio changes with the transmission rate r 
for different weight distributions. When links on a homogeneous network are distributed at 
random, the increase in the magnitude of one specific link weight (e.g. w±) accompanied 
by a decrease in its frequency leads to smaller R values. This is to be expected since the 
contribution of the different link types in this case is kept constant (piWi = P2W2 = 0.5) and 
this implies that the overall weight of the network links accumulates in a small number of 
highly weighted links with most links displaying small weights and thus making transmission 
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less likely. The statement above is more rigorously underpinned by the results of Theorem 
1 & 2 which clearly show that equal or more homogeneous weights lead to higher values of 
the basic reproductive ratio. For the case of fixed weight distribution, the changes in the 
value of i?o are investigated in terms of varying the weights, so that overall weight in the 
network remains constant. This is constrained by fixing values of p\ and P2 and, in this case, 
the highest values are obtained for higher values of w\. The flexibility here is reduced due 
to pi and P2 being fixed, and a different link breakdown may lead to different observations. 
The top continuous line in Fig. [TJ corresponds to the maximum Rq value achievable for both 
models if the piw± + P2W2 = 1 constraint is fulfilled. 

3.2 ii^-like threshold for the SIR model - a pairwise model per- 
spective 

To compute the value of Ro-like quantity from the pairwise model, we use the approach 
suggested by Keeling jH], which utilises the local spatial/network structure and correctly 
accounts for correlations between susceptible and infectious nodes early on in the epidemics. 
This can be achieved by looking at the early behaviour of [S7]i/[T] = Ai and [57] 2/ [i] = A2 
when considering links of only two different weights. In line with Eames [27], we start from 
the evolution equation of [I] 

[I] = (tw^SIMI] + tw 2 [SI] 2 /[I] - 7 ) [J], 

where from the growth rate twi\\ + TW2X2 — 7 it is easy to define the threshold quantity R 
as follows, 

™iAi + ™ 2 A 2 
R= . (11) 

7 

For the classic closure (J5J), one can compute the early quasi-equilibria for Ai and A 2 directly 
from the pairwise equations as follows 

Ai = and A 2 = ■ — . 

TW\ + 7.K TW2 + 7-Tl 

Substituting these into (TTTi) and solving for R yields 

D R1 + R2 + ^(Ri + R2) 2 + ±RJhQ nn , 
R = g ' ^ ' 

where 

rwi[(k - l)pi - 1] rw 2 [{k - l)p 2 - 1] 
Hi = , K 2 = , 

7 7 

fc-2 



Q 



p-i) Pl -i]p- i) P2 -iy 

with details of all calculations presented in Appendix D. We note that R > 1 will result in 
an epidemic, while R < 1 will lead to the extinction of the disease. It is straightforward 
to show that for equal weights, say W, the expression above reduces to R = rW(k — 2)/^ 
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which is in line with R value in [44J for unclustered, homogeneous networks. Under the 
assumption of a fixed total weight W, one can show that similarly to the network-based 
basic reproductive ratio, R achieves its maximum when W\ = w 2 = W. 

In a similar way, for the modified closure ([6]), we can use the same methodology to derive 
the threshold quantity as 



where 



D R1 + R2 + V( R i + R i) 2 + ^RiMQ - 1) 

K= n ) [ L6 ) 



TWi(k 1 -2) rw 2 {k 2 -2) k x k 2 
ti\ — , K 2 — , (4 



7 7 (h-2)(k 2 -2y 

For this closure once again, R > 1 results in an epidemic, while for R < 1, the disease dies 
out. Details of this calculations are shown in Appendix D. It is noteworthy that one can 
derive expressions ( lT2"j) and f[T3"j) by considering the leading eigenvalue of linearization of 
system flU) near its disease-free steady state with the corresponding pairwise closures given 
in (ED and ©. 

Finally, we note that this seemingly i?o-lookalike, R = rW(k — 2)/ r y for the equal weights 
case Wi = w 2 = W is a multiple of (k — 2) as opposed to [k — 1) as is the case for the 
Rq derived based on the individual-based perspective, where, for equal weights, R\ = R^ = 
rW(k— 1) / (rW+7). This obviously highlights the strong dependency of Rq on the modelling 
framework and also the difficulty in trying to reconcile findings based on different models. 

3.3 The performance of pairwise models and the impact of weight 
distributions on the dynamics of epidemics 

To evaluate the efficiency of the pairwise approximation models, we will now compare nu- 
merical solutions of models ([3]) and (J4]) to results obtained from the corresponding network 
simulation. The discussion around the comparison of the two models is interlinked with the 
discussion of the impact of different weight distributions/patterns on the overall epidemic 
dynamics. We being our numerical investigation by considering weight distributions with 
moderate heterogeneity. This is illustrated in Fig. EJ where excellent agreement between 
simulation and pairwise models is obtained. The agreement remains valid for both SIS 
and SIR dynamics, and networks with higher average link weight lead to higher prevalence 
levels at equilibrium for SIS and higher infectiousness peaks for SIR. 

Next, we explore the impact of weight distribution under the condition that the average 
weight remains constant (i.e. piW\ + p 2 w 2 = 1, where without loss of generality the average 
weight has been chosen to be equal to 1). First, we keep the proportion of edges of type 
one (i.e. with weight wi) fixed and change the weight itself by gradually increasing its 
magnitude. Due to the constraint on the average weight and the condition p 2 = 1 — pi, the 
other descriptors of the weight distribution follow. Fig. [3] shows that concentrating a large 
portion of the total weight on a few links leads to smaller epidemics, since the majority of 
links are low-weight and thus have a small potential to transmit the disease. This effect 
is exacerbated for the highest value of W\\ in this case 95% of the links are of weight 
w 2 = (1 —piWi)/(l —pi) = 0.5/0.95 leading to epidemics of smallest impact (Fig. Eta)) and 
smallest size of outbreak (Fig. EJ^b)). 
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While the previous setup kept the frequency of links constant while changing the weights, 
one can also investigate the impact of keeping at least one of the weights constant (e.g. the 
larger one) and changing its frequency. To ensure a fair comparison, here we also require 
that the average link weight over the whole network is kept constant. When such highly 
weighted links are rare, the system approaches the non-weighted network limit where the 
transmission rate is simply scaled by W2 (the most abundant link type). As Fig. H] shows, in 
this case, the agreement is excellent, and as the frequency of the highly weighted edges/links 
increases, disease transmission is less severe. 

Regarding the comparison of the pairwise and simulation models, we note that while the 
agreement is generally good for a large part of the disease and weight parameter space, the 
more extreme scenarios of weight distribution result in poorer agreement. This is illustrated 
in both Figs. [3] and H] (see bottom curves), with the worst agreement for the SIS dynamics. 
The insets in Fig. [3] show that increasing the average connectivity improves the agreement. 
However, the cause of disagreement is due to a more subtle effect driven also by the weight 
distribution. For example, in Fig. HJ the average degree in the network is 10, higher then 
used previously and equal to that in the insets from Fig. |3j but despite this, the agreement 
is still poor. 

The two different weighted network models are compared in Fig. This is done be using 
the same link weights and setting pi = k\/k and P2 = k-ijk. Epidemics on network with 
random weight distribution grow faster and, given the same time scales of the epidemic, this 
is line with results derived in Theorem 1 & 2 and findings concerning the growth rates. The 
difference is less marked for larger values of r where a significant proportion of the nodes 
becomes infected. 

In Fig. [H]the link weight composition is altered by decreasing the proportion of highly- 
weighted links. As expected, the reduced average link weight across the network leads to 
epidemics of smaller size while keeping the excellent agreement between simulation and 
pairwise model results. 

4 Discussion 

The present study has explored the impact of weight heterogeneity and highlighted that the 
added heterogeneity of link weights does not manifest itself in the same way as most other 
heterogeneities in epidemic models on networks. Usually, heterogeneities lead to an increase 
in Rq but potentially for final size to fall. However, for weighted networks the concentration 
of infectiousness on fewer target link, and thus target individuals, leads to a fall in Rq for 
both homogeneous random and fixed weight distribution models. Increased heterogeneity in 
weights accentuates the locality of contact and is taking the model further from the mass- 
action type models. Infection is concentrated along a smaller number of links, which results 
in wasted infectivity and lower Rq. This is in line with similar results [T5| [16], [70] where 
different modelling approaches have been used to capture epidemics on weighted networks. 

The models proposed in this paper are simple mechanistic models with basic weight 
distributions, but despite this they provide a good basis for analysing disease dynamics on 
weighted networks in a rigorous and systematic way. The modified pairwise models have 
performed well, and provide good approximation to direct simulation. As expected, the 
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agreement with simulations typically breaks down at or close to the threshold, but away from 
it, pairwise models provide a good counterpart or alternative to simulation. Disagreement 
only appears for extreme weight distributions, and we hypothesise that this is mainly due to 
the network becoming more modular with islands of nodes connected by links of low weight 
being bridged together by highly weighted links. A good analogy to this is provided by 
considering the case of a pairwise model on unweighted networks specified in terms of two 
network metrics, node number N and average number of links k. The validity of the pairwise 
model relies on the network being connected up at random, or according to the configuration 
model. This can be easily broken by creating two sub-networks of equal size both exhibiting 
the same average connectivity. Simulations on such type of networks will not agree with the 
pairwise model, and highlights that the network generating algorithm can push the network 
out of the set of 'acceptable' networks. We expect that this or similar argument can more 
precisely explain why the agreement breaks down for significant link-weight heterogeneity. 

The usefulness of pairwise models is illustrated in Fig.[7J where the I/N values are plotted 
for a range of r values and for different weight distributions. Here, the equilibrium value has 
been computed by finding the steady state directly from the ODEs ([3]) by finding numerically 
the steady state solution of a set on nonlinear equations (i.e. [A] = and [AB] = 0). To 
test the validity, the long term solution of the ODE is plotted along with results based on 
simulation. The agreement away from the threshold is excellent and illustrates clearly the 
impact of different weight distributions on the magnitude of the endemic threshold. 

The models proposed here can be extended in a number of different ways. One po- 
tential avenue for further research is the analysis of correlations between link weight and 
node degree. This direction has been explored but in the context of classic compartmental 
mean-field models based on node degree [121 [57]. Given that pairwise models extend to 
heterogeneous networks such avenues can be further explored to include different type of 
correlations or other network dependent weight distributions. While this is a viable direc- 
tion, it is expected that the extra complexity will make the pairwise models more difficult to 
analyse and disagreement between pairwise and simulation model more likely. Another the- 
oretically interesting and practically important aspect is the consideration of different types 
of time delays, representing latency or temporary immunity [13], and the analysis of their 
effects on the dynamics of epidemics on weighted networks. The methodology presented in 
this paper can be of wider relevance to studies of other natural phenomena where overlay 
networks provide effective description. Examples of such systems include the simultaneous 
spread of two different diseases in the same population [12] , the spread of the same disease 
but via different routes [UJ or the spread of epidemics concurrently with information about 
the disease [3"8| H5]. These areas offer other important avenues for further extensions. 
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5 Appendix 

5.1 Appendix A - Reducing the weighted pairwise models to the 
unweighted equivalents 

We start from the system 

[S]=l[I]-rEn=^n[SI] n , 

[SI] m = 7 ([//] m - [SIU + r£f =1 w n ([SSI] mn - [ISI] nm ) - rw m [SlU ( 14 ) 
[II] m = -2 7 [//] m + 2r£f =1 w n [ISI] nm + 2rw m [SI] m , 

[SS] m = 2j[SI] m - 2rY,n=i w n[SSI] mn , 

where m — 1,2, ... , M. To close this system of equations at the level of pairs, we use the 
approximations 

k-l [AB] m [BC] n 
[ABC] mn - — ^ . 

To reduce these equations to the standard pairwise model for unweighted networks we use 
the fact that J2m=i [AB] m = [AB] for A, B e {S, 1} and aim to derive the evolution equation 
for [AB]. Assuming that all weights are equal to some W, the following relations hold, 

M 

[si] = £ [si] m 

m=l 

M / M 

£ ^([ /J W - i SI U + rJ2M[SSI] mn - [ISI] nm ) - rw m [SI], 

m=l \ n=l 

M M 

= 7 ([//] - [SI]) - rW[SI] + rWj2J2 ([<^U " [ISI] nm ), 

m=l n=l 

where the summations of the triples can be resolved as follows, 

MM , , M M 



m=l n=l 



m=l n=l 



k-l[SS][SI] _ 

~k [s^~ [bbJl 



17 



Using the same argument for all other triples, the pairwise model for weighted networks 
with all weights being equal (i.e. W — 1) reduces to the classic pairwise model, that is 



[S]= 1 [I]-r[SI], 
[i]=r[5i]- 7 [i], 

E™ =1 [SI] = [SI] = 7 ([/i] - [SI]) + r[^7] - [ISI] - [SI], 

Tl=, [ir\ = [ii] = -Hm + MVsr\ + [si]), 

£™=i [55] = [55] = 2 7 [5il-2r[55il. 



A similar argument holds for the pairwise model on weighted networks with SIR dynamics. 
5.2 Appendix B - Proof of Theorem 1 

We illustrate the main steps needed to complete the proof of Theorem 1. This revolves 
around starting from the inequality itself and showing via a series of algebraic manipulations 
that it is equivalent to a simpler inequality that holds trivially. Upon using that p\k = ki, 
p 2 k = k 2 and p 2 + Pi = 1, the original inequality can be rearranged to give 



\AfM ~ !) r i - ( fc 2 - l)^] 2 + 4A;i/c 2 rir 2 < (h - l)n + (k 2 - l)r 2 + 2rip 2 + 2r 2 pi. (15) 



Based on the assumptions of the Theorem, the right-hand side is positive, and thus this 
inequality is equivalent to the one where both the left- and right-hand sides are squared. 
Combined with the fact that p 2 = 1 — p±, after a series of simplifications and factorizations 
this inequality can be recast as 



4pi(l - Pi){r\ + r\) + 8fcpi(l - Pi)nr 2 < Ak Pl (l - p x ){r\ + r\) + 8pi(l - Pi)nr 2 , (16) 



which holds trivially and thus completes the proof. We note that in the strictest mathemat- 
ical sense the condition of the Theorem should be (ki — l)ri + (k 2 — l)r 2 + 2r\p 2 + 2r 2 p\ > 0. 
This holds if the current assumptions are observed since these are stronger but follow from 
a practical reasoning whereby for the network with fixed weight distribution, a node should 
have at least one link with every possible weight type. 

5.3 Appendix C - Proof of Theorem 2 

First, we show that R^ is maximised when w\ = w 2 = W. Rq can be rewritten to give 



which can be further simplified to 



^(l-piXn-ra) 2 ^-!) >0, 



(17) 




(18) 
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Maximising this given the constraint W\p\ + w 2 (l — pi) = W can be achieved by considering 
i?Q as a function of the two weights and incorporating the constraint into it via the Lagrange 
multiplier method. Hence, we define a new function f(wi,w 2 , A) as follows 

f(wi,w 2 , X) — (k — 1) (pi^- + (I -Pi)- T "' 2 



' twi + r rw 2 + r 

+X(w 1 p 1 + w 2 (l - p x ) - W). 
Finding the extrema of this functions leads to a system of three equations 

df (k - l)piT7 



dwi (twi + 7) 2 
df (k-l)(l-p 1 )r 1 



+ Api = 0, 

+ A(l-pi) = 0, 



(19) 



dw 2 (tw 2 + 7) 

wipi + w 2 (l —pi) — W — Q. 
Expressing A from the first two equations and equating these two expressions yields 

{k — 1)t7 (k — l)r7 
(rwi+7) 2 (™ 2 + 7) 2 ' 
Therefore, 

Wl = w 2 = W, (20) 

and it is straightforward to confirm that this is a maximum. 

Performing the same analysis for Rq is possible but it is more tedious. Instead, we 
propose a more elegant argument to show that Rq under the constraint of constant total 
link weight achieves its maximum when w\ = w 2 = W. The argument starts by considering 
_Rq when Wi = w 2 = W. In this case, and using that r 2 = r 1 = r = rWj (rW + 7) we can 
write, 

2 , (h - l)n + (k 2 - l)r 2 + ^/[(h - l)n - (k 2 - l)r 2 ] 2 + ±k x k 2 r x r 2 

tin 



L 



r{h + k 2 -2) + ^[{h - I) - {k 2 - I)}* + 4r 2 hk 2 



2 

= r(2k 1 + 2k-2) =r(h + k2 _ l) = (k _ l)r 

However, it is known from Theorem 1 that Rq < Rq, and we have previously shown that Rq 
under the present constraint achieves its maximum when W\ = w 2 = W, and its maximum 
is equal to {k — l)r. All the above can be written as 

Rl<Rl<(k- l)r. (21) 
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Now taking into consideration that = (k — l)r, the inequality above can be written as 

R 2 < Rl < (k - l)r = Rl* , (22) 

and this concludes the proof. 

5.4 Appendix D - The i? -like threshold R 

Let us start from the evolution equation for [/](£), 

[I] = r(w 1 [SI} 1 + w 2 [SI} 2 )- 1 [I] 

= (tW 1 \ 1 + TW 2 X 2 ~ 

where Ai = ^jj 1 and A2 = ^jj^ an d let R be defined as 

R = TWlAl + TW2A2 . (23) 

7 

Following the method outlined by Keeling [H] and Eames [27], we calculate the early quasi- 
equilibrium values of Ai i2 as follows: 

\ 1 = o**[si\ l [i\ = [i\[si\ 1 , 

A 2 = ^[5/] 2 [/] = [/]][5/] 2 . 
Upon using the pairwise equations and the closure, consider [5/]i[7] = [I][5i]i: 
[SIUI] = (rw^SSIln + rw 2 [SSI] 12 - T Wl iISI] n - rw 2 [ISI] 21 - rw 1 [SI] 1 - j[SI]x)[I] 

= (Tw 1 [Sl\ 1 + TW 2 [SI\ 2 -'y[I\)[SI\ 1 . (24) 
Using the classical closure 

k — 1 [AB]i[-BC] 2 



[ABC] 12 
[ABC} 21 



k [B] 

k-ljABUBC^ 
k [B] 



and making the substitution : [SI] 1 = Xi[I] , [SI] 2 = X 2 [I], [I] < 1, [5] « N, [55] 1 w fciVpi, 
[55] 2 ~ fciV(l — Pi) together with 7/? = tu^Ai + rw; 2 A 2 , we have 

(rwxXi + rw 2 X 2 )kp 1 - (rwiAi + ™ 2 A 2 )pi - (tiuiAi + ™ 2 A 2 )A X - riOiAi = 0, 

which can be solved for Ai to give 

= j(k - l) Pl R 

TW\ + ^R 



20 



Similarly, A 2 can be found as 



rw 2 + 7-R v ; 



Substituting the expressions for Ai j2 into the original equation for R yields 

D A + B + + Bf + AT 2 wiw 2 {k - 2) 

2 7 

where A = TW\[(k — l)p\ — 1] and B = rw 2 [(fc — l)p 2 — 1]. If we define 

= rwi[(fc - - 1] ^ d = Tw 2 [{k - l)p 2 - 1] 
1 7 2 7 

the expression simplifies to 

7?i + 7? 2 + ^(#1 + R2) 2 + ^R 2 Q 



where Q 



R = 
(k-2) 



[(k — l)pi - 1] [(k — l)p 2 — 1] 
Substituting the modified closure 

[ABC] U = 

[ABC] 12 = 

[ABC] 21 = 



[ABC] 
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fci-l[AB]![BC]i 
fei [5] : 

[AB]![EC] 2 
[B] ' 

[ABjajgCJi 
[S] ' 

k 2 - 1 [AB] 2 [BC] 2 
fc 2 [B] 



into 021]) and making further substitution : [SI]i = Ai[I], [57] 2 = A 2 [7], [I] < 1, [5] « iV, 
[SS]x « fciJV, [55] 2 ~ feTV, we have 

(riuiAi + rw 2 \ 2 )ki — (tiuiAi + rw 2 A 2 )Ai — 2TU7A1 = => Ai ' 



2rw! + 7i2 

Similarly, the equation [57] 2 [7] = [7] [57] 2 yields 
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Substituting these expressions for Ai j2 into f )23|) . we have 
r(wiki + w 2 k 2 ) - 2t(w 1 + w 2 ) 



R 
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[2t(wi + w 2 ) - r(w 1 k 1 + w 2 k 2 )] + 8r 2 WiW 2 (fci + ^2 - 2) 

M 5f 

If we define 

rwi(A;i - 2) ru^fe - 2) 

itl — , K 2 — , 

7 7 

the above expression for R simplifies to 



D R1 + R2 + ^(Ri + R 2 ) 2 + ^RiR 2 (Q-l) , ntt . 
R = (2o) 



where 

Q hh 



(k 1 -2)(k 2 -2)' 
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Figure 1: Plots of the basic reproductive ratio R for the two types of weighted networks 
with different weight and weight frequency combinations, but with p±Wi + P2W2 = 1. The 
case of homogenous networks with weights assigned at random (continuous lines with p\ and 
Wi given in the legend) considers the situation where the contribution of the two different 
weight types is equal (piu>i = P2W2 = 0.5) but with weight W\ increasing (top to bottom) and 
its frequency decreasing. Increasing the magnitude of weights but reducing their frequency 
leads to smaller Rq values. The case of homogeneous networks with fixed number of links 
of type W\ and W2 (dashed lines with only w\ given in the legend) illustrates the situation 
where W\ increases (bottom to top) while pi = k\/k = 1/3 andp2 = (k—ki)/k = 2/3 remain 
fixed. Here the opposite tendency is observed with increasing weights leading to higher R 
values. Finally, for the randomly distributed weights case, setting p 1 = 1/3, W\ = 1.4 and 
observing p±Wi + P2W2 = 1, we obtain _R (*) values which compare almost directly to the 
fixed-weights case (top, dashed line). Other parameters are set to k = 6, k± = 2 and 7 = 1. 
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Figure 2: The infection prevalence (I/N) from the pairwise and simulation models for ho- 
mogeneous random networks with random weight distribution (ODE: solid line, simulation: 
dashed line and (o)). All nodes have degree k = 5 with N = 1000, Io = 0.05iV, 7 = 1 and 
r = 1. From top to bottom, the parameter values are: w\ = 5,pi = 0.2, w 2 = 1.25, p 2 = 0.8 
(top), and Wi = 0.5, pi = 0.5, w 2 = 1.5, p% = 0.5 (bottom). The left and right panels 
represent the SIS and SIR dynamics, respectively. 
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Figure 3: The infection prevalence {I/N) from the pairwise and simulation models for ho- 
mogenous networks with random weight distribution (ODE: solid line, simulation: dashed 
line and (o)). All numerical tests use N = 1000, Iq = 0.05iV, k — 5, 7 = l,r = 1 
and pi = 0.05 (p 2 = 1 — Pi = 0.95). From top to bottom, w\ = 2.5,5,10, w 2 = 
0.875/0.95,0.75/0.95,0.5/0.95. The weight distributions are chosen such that the average 
link weight, p\Wi + p 2 w 2 = 1, remains constant. Insets of (a) and (b): the same parameter 
values as for the lowest prevalence plots but, with k — 10 and r = 0.5. The left and right 
panel represent the SIS and SIR dynamics, respectively. 
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Figure 4: The infection prevalence {I/N) from the pairwise and simulation model for ho- 
mogenous networks with random weight distribution (ODE: solid line, simulation: dashed 
line and (o)). All numerical tests use N = 1000, Iq = 0.05iV, k — 10, 7 = l,r = 0.5 and 
Wi = 10. From top to bottom, P{w x ) = 0.01,0.05,0.09, w 2 = 0.9/0.99,0.5/0.95,0.1/0.91. 
Here also p 2 = 1 — p\ and p\W\ + p 2 w 2 = 1. The left and right panel represent the SIS and 
SIR dynamics, respectively. 
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Figure 5: The infection prevalence (I/N) based on random (model 1) and fixed (model 2) 
weight distribution (ODE: black (1) and blue (2) solid line, simulation results: same as ODE 
but dashed lines, and (o) and (*)). All numerical tests use N = 1000, Iq = 0.05iV, k = 10, 
k\ = 2, k 2 = 8, pi = ki/k, p 2 = k 2 /k, w\ = 10, w 2 = 1.25 and 7 = 1. The rate of infection 
r = 0.5 (top) and r = 0.1 (bottom). The left and right panel represent the SIS and SIR 
dynamics, respectively. 
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Figure 6: The infection prevalence {I/N) for a fixed weight distribution (ODE: solid lines, 
simulation results: dashed lines and (o)). All numerical tests use N = 1000, Iq = 0.05iV, 
k — 6, 7 = 1, r — 1 and w± = 1.4, w 2 = 0.8. From top to bottom : k\ = 5,4,3,2,1 and 
k% — k — k\. The left and right panel represent the SIS and SIR dynamics, respectively. 
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Figure 7: Endemic steady state from the SIS model on networks with random weight 
distribution. The continuous lines correspond to the steady state computed numerically by 
setting all evolution equations in the pairwise system to zero. These are complemented by 
finding the endemic steady state through direct integration of the ODE system for a long- 
enough time (o), as well as direct simulation (*). The first marker corresponds to r = 0.3 
followed by r = 0.5, 1.0, . . . , 3.0. All results are based on: k — 5, 7 = 1 and W\ — 10, Wi — 1. 
From top to bottom : p± = 0.9, 0.5, 0.1, 0.01 and p% — 1 — Pi- 



